HEAVY MESON DYNAMICS 
IN A QCD RELATIVISTIC POTENTIAL MODEL 



FULVIA DE FAZIO 

Dipartimento di Fisica dell'Universitd di Bari, 
Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 
Via Amendola 173, 70126, Bari, Italy 

ON 

ON ■ ABSTRACT 

We use a QCD relativistic potential model to compute the strong coupling constant 
D ' g appearing in the effective Lagrangian which describes the interaction of — and 

C/3 , 1~ qQ states with soft pions in the limit niQ — > oo. We compare our results with 

' other approaches; in particular, in the non relativistic limit, we are able to reproduce 

. the constituent quark model result: g = 1, while the inclusion of relativistic effects 

due to the light quark gives g — g, in agreement with QCD sum rules. We also 
estimate heavy meson radiative decay rates, with results in agreement with available 
experimental data. 

Q\ \ 1. The Strong Coupling Constant gp,*D-K 
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The decay D* + — > D°tt + is described in terms of a strong coupling constant 
9d*d-k defined by: < D°(k)n + (q)\D* + (p, e) > = go*Dn CLEO collaboration 

measurement a: BR(D* + — > D°n + ) ~ 68.1 ± 1.0 ± 1.3% and the upper boundcl : 
(JL, ■ T(D* + ) < 131 KeV provide us with the constraint : go*D-K < 20.6. 

The interest in the evaluation of this coupling constant is manifold. The form 
factor Fi(q 2 ) describing the semileptonic decay B — > tt£u is believed to be dominated 
by the B* pole, so that its value at maximum transferred momentum is proportional 
to gB*Bir, which is related to gD*Dn by i: gp*p n = 2 nip g/f^ , where P(P*) is a 0~ 
(1~) heavy meson of mass nip, and g is independent of nip. Besides, g appears in 
the effective Lagrangian describing the interaction between heavy mesons and light 
Nambu-Goldstone bosons 

Big 

In non relativistic quark models ill g ~ 1, while recent QCD sum rules I and 
HQET §0 analyses give: g ~ 0.2 — 0.4. We wish to show that the inclusion of 
the relativistic effects in the bound state can lower the value g — 1, providing an 
explanation of the discrepancy between the different approaches. 

We shall obtain g in the framework of a QCD relativistic potential model In 
this model the qQ heavy states D and D* are described in terms of the creation 
operators of the constituent quarks and of a meson wave function ip, normalized 
according to: -^p- / dk\tp\ 2 = 2^Jm 2 D + p 2 , where p is the meson momentum and 
k is the quark relative momentum, ip satisfies a Salpeter equation* which includes 



*This equation arises from the bound-state Bethe Salpeter equation by considering the instantaneous 
approximation and restricting the Fock space to the qQ pairs 113. 



relativistic effects in the quark kinematics. In this equation the interquark potential 
V coincides, in the meson rest frame, with the Richardson potential, with a linear 
behaviour for large distances, reproducing QCD confinament, and a coulombic shape 
for small distances, with the further assumption that V is constant near the origin in 
order to avoid unphysical singularities. The values of the quark masses, obtained by 
fits to meson masses, are: m u = rrid = 38 MeV, m s = 115 MeV, m c = 1452 MeV, 
m b = 4890 MeV. 

To evaluate go*D-K let us consider the matrix element of the axial current 
between the states D* and D: 

<D°(k)\A ll \D* + (p,e)>=-t{ e M K+m )4((| 2 ) 

- — 2m D ,q^\Az(q 2 ) 

where 2mD*A 3 = {mo + mr>*)Ai + (mo* — rriD)A<i (q = p — k). 

Taking the derivative of A^, we can link the l.h.s. of Eq. (|]) to the matrix 
element of the pseudoscalar current between the same states, which is supposed to 
be dominated by the 7r + pole. Performing the overlap of the states D, D* and the 
axial current, we have in the chiral limit: gD*Dw — 2m f D * A (0), and finally: 

1 r 00 

g = A (0) = - / dku(k) 
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where E q = ^Jk 2 + m 2 q and u(k) = ■ The wave functions u{k) come from a 
numerical solution of the Salpeter equation. We obtain 0: A (0) = 0.4 (D case); 
A (0) = 0.39 (B case), giving so: Qd*d-k = 12.3 and gs-B-K = 31.7. We can observe 
that we have only 2% deviation from the scaling result gD*Dn/gB*Bn = rrio/mB- 

It is interesting to notice that, in the non relativistic limit, i.e. E q ~ m q ^> k, we 
obtain: g = /o°° dk\u{k)\ 2 = 1, reproducing the constituent quark model result. 
On the other hand, in the limit m q — > 0, mg — > oo, the result is: g = 1/3, in 
agreement with the QCD sum rules determination. 
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2. Radiative Heavy Meson Decays 

The evaluation of radiative decay rates involves the knowledge of the matrix ele- 
ment of the electromagnetic current J^ m ' between the states D* and D: 

< D+(k)\j;- m -\D*+(p,e) >= hi + ^) e, vaP e^V , (3) 

where cq (e q ) is the heavy (light) quark electric charge. In the framework of the 
relativistic QCD model we are able to evaluate the "effective" masses A q , Aq by 
computing the overlap of meson states and the electromagnetic current. We find 



Table 1. Heavy meson radiative decay rates 



Decay rate/ BR 


theory 


experiment 


r(£)*+) 




46.21 ifeV 


< 131 KeV 


BR(D* + 


-> D + 7T ) 


31.3% 


30.8 ±0.4 ±0.8% 


BR(D* + 


-> D°7T+) 


67.7% 


68.1 ±1.0 ±1.3% 


BR{D* + 




1.0% 


1.1 ± 1.4 ± 1.6% 


T(D*°) 




41.6 KeV 




BR(D*° 


-> £)°7T ) 


50.0% 


63.6 ± 2.3 ± 3.3% 


BR(D*° 


^/J 7 ) 


50.0% 


36.4 ±2.3 ±3.3% 


T(D* S ) = 


r(L>: - d s1 ) 


0.382 KeV 




T(B*-) = 


= T(B*~ -> 5-7) 


0.243 KeV 




T(B*°) = 


= 1(5*° -> 5° 7 ) 


9.2 10~ 2 fTeF 




T(B* S ) = 


r(5; - £? s7 ) 


8.0 10~ 2 KeV 





A c ~ m c = 1.57 GeV, A. ~ m 6 = 4.95 GeV and A q = 0.48 GeV > m q . The B.R. 
are reported in Table I cB together with the available experimental data. 

We conclude that the inclusion of relativistic effects in a QCD potential model has 
allowed us to explain the discrepancy among different models in the evaluation of the 
heavy meson coupling constant with soft pions. Moreover, the same model, applied 
to heavy meson radiative transitions, gives results in agreement with experimental 
data. 

I thank P. Colangelo and G. Nardulli for their precious collaboration. 
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